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ABSTRACT 



We modified the Press-Schechter (PS) formalism and then analytically 
derived a constrained mass distribution function n(M\ip) for the regions having 
some specified value of the primordial gravitational potential, <p. The resulting 
modified PS theory predicts that gravitationally bound clumps with masses 
corresponding to rich clusters are significantly biased toward the regions of 
negative primordial potential - the troughs of the potential. The prediction is 
quantitative, depending on the mass and the depth of the troughs, which can 
be tested in large N-body simulations. As an illustration of the magnitude of 
the effect we calculate the constrained mass function for the CDM model with 
T = Qh = 0.25 normalized to as = 1. In particular, we show that the probability 
of finding a clump of mass 10 14 — 10 15 /i~ 1 M Q in the region of negative initial 
potential is 1.3 — 3 times greater (depending on the mass) than that in the 
region of positive initial potential. The scale of the potential fluctuations 
Rip = y/Scr^/a^ is shown to be « 120/i _1 Mpc for the spectrum in question. The 
rms mass density contrast on this scale is only about <ys{R ip ) ~ 0.03. Assuming 
that the modified PS theory is statistically correct, we conclude that clusters are 
significantly biased {b > 10, b is a bias factor defined by An c //n c / = bAp m /p m ) 
toward the regions having negative initial potential. 

Subject headings: cosmology: theory — large-scale structure of 
universe 
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1. INTRODUCTION 

Assuming the standard hierarchical model of the structure formation from Gaussian 
fluctuations due to gravitational instability, we study the effect of primordial gravitational 
potential fluctuations on massive objects such as galaxy clusters and perhaps superclusters, 
i.e., clusters of clusters (Pahcall fc Soneira 1984Q . We employ the PS formalism as a tool 
and modify it for this study. 

Some effect of the primordial gravitational potential upon the structure formation 
has been already noted. [Kofman fc Shandarin (1988)| have noticed that the adhesion 
approximation predicts that the formation of voids is associated with positive peaks 
of the primordial gravitational potential. Sahni, Sathyaprakash, & Shandarin (1994) 
studied the effect and measured a significant correlation between the sizes of voids and 
the value of primordial gravitational potential in numerical simulations of the adhesion 
model. By investigating the evolution of correlation between the potential and the density 
perturbations, Buryak, Demianski, & Doroshkevich (1992) showed that the formation 
of super large scale structures is mainly determined by the spatial distribution of the 
gravitational potential. Recently, Madsen et al. (1997) have demonstrated by N-body 
simulations that the under dense and the over dense regions are closely linked to the regions 
with the positive and the negative gravitational potential respectively. Thus, given all these 
results showing the important role of the primordial gravitational potential in the structure 
formation, it would be interesting to calculate the effect of the primordial potential upon 
the mass distribution function. 

The mass distribution function n(M) is defined such that n(M)dM is the comoving 
number density of gravitationally bound objects in the mass range (M,M + dM). The 
standard Press-Schechter (hereafter, PS) formalism provides an effective tool to evaluate 
n(M) in spite of various criticism on it (see [Monaco 1998| ), and is widely used in cosmology 
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(e.g. Gross et al. 1997 ; Kitayama fc Suto 1997 ; Bahcall fc .ban 199S ; Robinson, Gawiser, & 
|Silk 1998| ; [Wang fc Steinhardt 1998| ). Also, |Lee fc Shandarin ( 1998 )| have shown by applying 
the dynamics based on the Zel'dovich approximation to the PS formalism that it is very 
robust with respect to the underlying dynamics. 



The following two equations represent the essence of the PS formalism (Press fc 

dF 



Schechter 1974|) : 



n(M) = £ 



F(M) 



dM 1 
p(5)d5. 



(1) 
(2) 



Here p(5) is the probability density distribution of the linearly extrapolated density 
contrast 5 smoothed on a comoving filtering scale R which is related to the mass by 
M = M(R) = apR? . The proportionality constant a is either determined by the shape of 
the smoothing window function or sometimes is used as a free parameter in order to provide 
a better agreement with numerical results. In the case of a sharp k-space filter which is 
actually consistent with the PS formalism (see Peacock fc Heavens 1990|) , the filtering scale 
k c = 2n/R in k-space and mass are related as M = 67r 2 pk~ 3 . The density threshold value 
S c for collapse was originally given as 5 C ~ 1.69 according to the Top Hat spherical model. 
However, it has been shown that the lowered value of S c in the range from 1.3 — 1.6 gives a 
better fit in N-body simulations, which depends on the the initial spectrum and the type of 
the filter (e.g., |Gross et al. 1997|). 



In this Letter we investigate and show how much the primordial gravitational potential 
ip affects the mass distribution function of galaxy cluster. Modifying the PS formalism, we 
derive a constrained mass distribution function n(M\tp) defined as the comoving number 
densities of clumps of mass M in the regions where the primordial gravitational potential 
fluctuation satisfies some specified conditions. The Cold Dark Matter model (CDM) with 
T = Qh = 0.25 and a$ = 1 is used to demonstrate the significance of the effect. 
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2. 



MODIFICATION OF THE PS FORMALISM 



In order to incorporate the primordial gravitational potential fluctuations term 
into the above equations, we first derive the conditional probability density distribution 

p(S\(f < -tp c ) (y? c > 0): 



Here erf, er^, and a 2 are the density, velocity and the potential variances respectively; 
k =< 5 ■ <p > /(T$a v = <Tv/< 7 s< 7 <p 1S the crosscorrelation coefficient of the the density contrast 
S smoothed on the scale k c and the primordial (unsmoothed) potential fluctuations ip. The 
probability density distribution function for ip > <p c has the same form as equation (3) 
except for the opposite sign in front of the second error function term. 

In the case of a sharp k-space filter assumed throughout this Letter, erf, er^ and a 2 are 
given as 



where P(k) is the density power spectrum. It is worth stressing that the variances erf and a 2 
depend on k c , but the variance of the primordial potential, a 2 does not. We are interested 
in evaluating the effect of the primordial potential fluctuations and thus take er^ with no 
filtering. The long wave cutoff k\ xs 3000/i _1 Mpc (eq. 6) corresponds to the assumption 
that the waves longer than the cosmological horizon are irrelevant to the processes on scales 
of the structures in the universe (i.e. smaller than a few hundred of /i _1 Mpc). A particular 




(3) 




(4) 



(6) 



(5) 
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choice of k\ is unimportant since for the Harrison-Zel'dovich spectrum (assumed here) a 2 ^ 
depends on l/k t only logarithmically. 

As a result of the incorporation of the potential into the probability density distribution, 
the volume fraction F(M\(p < —(p c ) is now a function of both as and a v each of which in 
turn is a function of mass M: 

pod 

F{M\tp < -<p e ) = F[a s (M), a v (M)\ = / d8 p(5\<p < -<p c ). (7) 
Thus, equation (1) for the conditional mass function n(M\ip < —(p c ) becomes 



n(M|v? < -<p c ) = ^ 



dF das dF da v 



das dM da v dM 



For the initial power spectrum P(k), we use the fit given by Pardeen et al. (1986 



with T = Qh = 0.25 and a$ = 1 normalization. This choice of the parameters is in general 
agreement with both the COBE measurements and the galaxy two-point correlation function 
( [Peacock fc Dodds 1994[ ). Calculating das tV /dM numerically, we evaluate n(M\(p < — ip c ) 
through equations (4) to (8). We reserve the detailed description of the calculation for a 
companion paper in preparation. 

Fig. 1 shows the magnitude of the effect in terms of the cumulative mass functions. 
The upper panel shows the cumulative conditional mass functions n(> M\(p < —f c ) an d 
n(> M\tp > <f c ) (for each case of <p c = a v , 0) along with the unconditional PS mass function, 
n(> M). The cumulative mass functions were obtained by integrating the mass functions 
n(M\(p < — (p c ) numerically. The shaded area shows the fit to the observational mass 
function of rich clusters given by Bahcall fc Cen (1993) . The ratios of the conditional mass 



functions to the unconditional one are plotted in the lower panel of Fig. 1. 

We also calculate the probability that a clump with mass M is located in the potential 
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regions satisfying the chosen condition, for instance, ip < —(p c 

p(, e < -, e m = n(M %;' Pc) p(* < -9.), (9) 

where P(ip < —f c ) is the fraction of space satisfying the given condition (see Fig. 2). The 
other probabilities P(—<p c < V < 0|M), P(0 < <p < <p c \M), and P(y? c < <fi\M) can be 
obtained in a similar manner. 



3. RESULTS AND DISCUSSION 

The PS formalism has been proved to be a simple but very effective tool widely used 
for constraining cosmological models. We have modified it by considering the dependence of 
mass function on the initial perturbation of gravitational potential. The resulting modified 
PS theory predicts that the clumps with masses greater than roughly 10 1A h~ l M & have 
a noticeable tendency to form in the troughs of the primordial gravitational potential 
(the regions where the primordial potential fluctuations were negative). This quantitative 
prediction can be tested in large N-body simulations. Regardless of the outcome it will shed 
light on the PS formalism; if our prediction is confirmed, it will show a new potency of the 
PS technique. Otherwise a new limitation to the formalism will be established. 

Assuming that the prediction is correct at least qualitatively, [] we would like to discuss 
some of its obvious consequences. The scale of the initial potential 



f°°dkk~ 2 P(k) 

R v = vW<V = \ 3 joo dkp ^ « UOh^Mpc (10) 

does not depend on any ad hoc scale; the dependence on k\ is exremely weak (oc yjln{l/ki) 
for the Harrison-Zel'dovich spectra assumed here). It is, perhaps, worth mentioning that 



1 N-body simulations (e.g., Madsen et al. 1997) and the adhesion model (Sahni et al. 
1994) have already visually demonstrated this bias effect of the gravitational potential. 
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the scale of the potential is also practically independent of the smoothing scale unless it 
exceeds the value of a few tens of /i _1 Mpc. The density scale R& kc is determined by the 
scale of the smoothing window function k c that has only one "natural" scale corresponding 
nonlinearity k c = k n i. For the model in question the scale of the primordial potential 
is found to be R v ~ 120/i _1 Mpc. The scale of the density contrast field reaches this 
value Rs = ySag/as' ~ 120/i -1 Mpc only after it is smoothed on k c ~ 0.017/iMpc -1 . The 
corresponding density variance on this scale is ^(O.OlT/iMpc -1 ) « 0.03. On the other hand, 
the number of clumps with masses 10 14 — 1O 15 /i _1 M can easily be 30% greater in the troughs 
of the potential than the mean density n(> M) = 0.5[n(> M\tp < 0) + n(> M\tp > 0)] (see 
Fig. 1). Thus, the bias factor b (defined by the relation An cl /n cl = bAp m / p m ) reaches at 
least 10 on the scale about 120/i _1 Mpc. 

Qualitatively the bias phenomenon can be explained as follows. The initial density 
contrast is proportional to the Laplacian of the initial potential (5 oc V 2 <^). Therefore 
the two fields are cross-correlated: the positive peaks of 5 are more likely to be found 
in the troughs of the potential where it is negative. The correlation is not very 
strong (for k c = 0.25/iMpc -1 corresponding to a$ = 1 the crosscorrelation coefficient 
hi = o" 2 (0.25/iMpc _1 )/cr ¥ ,cr,5(0.25/iMpc _1 ) w 0.12). But the clusters are extreme objects 
corresponding to the tail of the mass function, and thus very sensitive to the environment. 
That is why the clusters put one of the strongest constraints on cosmological models 
( IKlypin fc Rhee 1994| , [Bond fc Myers 1996| , |Fan, Bahcall, fc Cen 19971 , |Bahcall fc Fan 1998| ). 



Incorporating the motion of mass into dynamics can only increase the bias effect due 
to the nonlinear effects although they are quite small on the scale in question. But, the 
point is not in the magnitude of the nonlinear effects but rather in their sign. On the scale 
of the potential the mass moves from the peaks of the potential to the troughs. Using the 
Zel'dovich approximation one can easily estimate the rms displacement of the mass on the 
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scale of the potential ( [Shandarin 1993Q : 



v ; ~3A _1 Mpc. (11) 



It is relatively small compared to the scale of the potential but coherent on the scale of 
the potential field, and therefore it can only enhance the bias effect. Another nonlinear 
effect is related to the rate of growth of perturbations. For the perturbations on the scale 
of a few Mpc the potential troughs/peaks may be viewed as patches with slower/faster 
expansion rate that corresponds to the increase/decrease of the rate of growth of small-scale 
perturbations. Similarly, the bias is enhanced in the redshift space because the velocity 
field is directed toward the troughs and away from peaks of the potential. Both effects can 
increase the bias by about 5% depending on the initial spectrum. 

Another way of calculating the constrained mass function would be using the 
peak-background split technique suggested by [Kaiser (1984)| to explain the enhanced 
correlation function of reach clusters. Obviously, the initial potential resembles the 
smoothed initial density field if the filter has a sufficiently large scale, but the former is 
never identical to the latter. The potential itself can be viewed as a smoothed density 
field with a very soft scale-free filter W(k) oc k~ 2 . Typically the density field is filtered 
with much harder filters (e.g. top-hat, Gaussian, or sharp fc-space filters), that impose 
the scale which is an ad hoc parameter. The magnitude of the bias in our approach is 
determined by the crosscorrelation of the density contrast smoothed at the scale (k c ) of 
nonlinearity {crs kc = 1) with the initial potential that does not have any ad hoc parameters. 
Probably, the value of the crosscorrelation coefficient determines the bias in the peak-split 
approach as well. The crosscorrelation of the density field 5k c smoothed on the scale of 
nonlinearity (k c = 0.25/iMpc^ 1 ) with the field 5^ smoothed on the scale of the potential 
(hp = 0.017/iMpc -1 ) is about 4 times weaker than the correlation of 5k c with the initial 
potential cp. Thus, we expect that the bias of galaxy clusters on such large scales as the 
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scale of the initial potential (~ 120/i _1 Mpc) is stronger toward the troughs of the potential 
fluctuations rather than to the peaks of the density fluctuations 5k c smoothed with the 
corresponding filter. We have not applied the split peak-background approach because 
it is not clear how to avoid arbitrarines in choosing the scale that splits the density into 
small-scale peaks and large-scale background field. This question requires a separate study. 

Applying this effect to observations one has to take into account the following issues. 
The gravitational potential does not evolve much on large scales especially in the Einstein-de 
Sitter universe (Kofman & Shandarin 1988; |Pauls fc Melott 1995| ; [Melott, Sathyaprakash 



fc Salmi 1996| ). Therefore, the potential at present is very similar to the primordial one on 



scales much greater than the scale of nonlinearity. A simple explanation to this in the frame 
of the standard scenario of the structure formation is due to the fact that the mass has 
been displaced by the distance about 10/i _1 Mpc (|Shandarin 1993| ). Therefore, the potential 
on scales greater than, say, 30/i _1 Mpc has been changed very little. 

Clusters can be used as statistical tracers of the potential. In addressing this question 
it is worth noting that the shot noise is an important factor since clusters are rare objects. 
Using the observational mass function (Bahcall & Cen 1993) one can estimate that an 
average spherical patch of the radius ~ 60/i -1 Mpc contains about 30 clusters with the 
masses greater than 10 14 h~ l M & . Thus, the shot noise is about 18% on this scale which is 
comparable with the bias itself (see Fig. 1, the bottom panel). However, Fig. 2 suggests 
that the most massive clusters [M > 1O 15 /i _1 M ] are very likely to reside in the regions of 
negative potential (P > 75%) and very unlikely in the regions of high potential (P < 5% if 
ip > a^). More detailed analysis will be present elsewhere. 

Probably, the best candidates for the markers of the troughs in the field of the 
primordial potential fluctuations are superclusters (defined as clusters of clusters) ( [Bahcall 



|fc Soneira 1984| ) especially with highest density enhancements (Shapley supercluster) and 



the giant geometrical patterns in the cluster distribution ( 



l'ully et al. 19921) . 
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Fig. 1. — In the upper panel the conditional cumulative mass function satisfying chosen 
potential condition is plotted. The solid, the long dashed, the dot-dashed, and the dashed 
lines correspond to the conditions ip < — o^, ip < 0, <p > 0, and <p > respectively, while the 
dotted line represents the unconditional cumulative PS mass function. The shaded area is la 
fit to the observational cumulative mass function of rich clusters by Bahcall and Cen (1993). 
In the lower panel the ratio of the conditional cumulative mass functions to the unconditional 
one is plotted for each condition. The CDM spectrum with T = 0.25 normalized to <r 8 = 1 
has been used. 
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Fig. 2. — The probability that a clump with mass M can be found in the regions satisfying 
chosen potential condition is plotted. The heavy solid, the heavy dashed, the solid, the 
dashed, the long dashed, and the dot-dashed lines correspond to the condition y? < 0, <p > 0, 
ip < —a v , —a v < ip < 0, < ip < a v , and ip > respectively. 



